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I. INTRODUCTION 

Curved space-time forms the basis for most classical theories of gravity, such as general relativity. These theories 
are usually based on a metric for four dimensional space-time jl| . Some of the basic concepts used in general relativity 
and related theories are transformation rules for tensors, the affine connection, and the relation of the metric to proper 
time along an observer's world line. A useful, but little-used concept, is that of the world function of space-time, as 
developed by J. L. Synge. The world function is essentially one-half of the squared measure between two points in 
space-time. The utility of the world function comes from the fact that it is closely related to experiments, and that it is 
a type of scalar quantity. Since the world function transforms as a kind of scalar, it allows us to formulate geometric 
quantities in a covariant way. Hence, the world function is a valuable tool for understanding the geometric ideas 
in metric theories of gravity, in three dimensional differential geometry and tensor analysis, and wherever arbitrary 
coordinate systems are used. As an example of the utility of the world function, I present its application to the 
problem of navigation in a curved space-time. This application actually goes beyond a simple pedagogical example 
because it deals with the real-world need for precise navigation and time dissemination. 

Consider the problem of an observer who wants to navigate in a curved space-time with respect to electromagnetic 
beacons. I use the word navigate to mean that the observer determines his or her coordinate position and coordinate 
time along their world line, in some system of space-time coordinates. I assume that the electromagnetic beacons 
continuously broadcast their space-time coordinates and that this information is imbedded in the emitted electro- 
magnetic signals. Furthermore, I assume that an observer at unknown coordinates (i ,x ) simultaneously receives 
these signals from four beacons. The observer's navigation problem is to compute his position (£ ,x ) from the four 
received emission event coordinates (t s ,x s ), s =1,2,3,4, of the electromagnetic beacons. 

In the case of fiat space-time, the observer must solve the four simultaneous equations 

|x -x s | 2 -c 2 (t -t s ) 2 = 0, s = l, ...4 (1) 

Equation (Q) contains four invariant statements: light signals travel 'on the light cone' from each emission event 
to the observer, (see Figure [l]). To resolve the branches of the light cones, the causality conditions t Q > t s , for 
s =1,. . .4, must be added. The relations in Eq. (|l|) are the basic navigation equations applied by users of satellite 
navigation systems, such as the U.S. Global Positioning System (GPS) and the Russion Global Navigation Satellite 
System (GLONASS) The coordinates of the four events, (i s ,x s ), correspond to particular radio emissions by 

the satellites. The emission event coordinates can be extracted from information transmitted by digital codes. 

Equation (|l|) is commonly used in two different ways. First, an observer may receive radio signals from satellites and 
compute his or her space-time coordinates (t a , x ) in terms of four known satellite emission events (t s , x s ). The second 
use of Eq. ([!]) is to locate a satellite, which is at an unknown position (i ,x ) in terms of four ground observations 
at known coordinates (t s ,x s ). For this use, we apply the causality condition t s > t Q , s =1,. . .4. In Eq. ( |l|), the 
assumption is made that space-time is flat and the speed of light c is constant. Furthermore, by using Eq. (uh we 
make the geometric optics approximation that the wavelength of the electromagnetic waves is small compared to all 
physical dimensions of the receiver and transmitter systems |||| . 

In recent years, there have been significant improvements in the stability of frequency standards and measurement 
techniques (t]||- Consequently, over satellite-to-ground distances, precise measurements should be interpreted within 
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the framework of a curved space-time theory [p|-|l2|] . Furthermore, the equations for navigation in space-time should 
be manifestly covariant and also invariant [|13| . 

In this work, I write down a generalization of the navigation Eq. ([!]) for curved space-time and give the detailed 
equations that must be solved for navigating in the vicinity of the Earth, both in Schwarzschild coordinates and 
in rotating coordinates. I still retain the geometric optics approximation, however, I take into account deviations 
from flatness to first order in the metric. This means that the flat-space light cones in Eq. (0) are replaced by 
equations for null geodesies. The required navigation equations are simply expressible in terms of the world function 
developed by J. L. Synge juj. The resulting formalism takes into account the delay of electromagnetic signals due 
to the presence of a gravitational field. The detailed equations have application to a user who wants to accurately 
compute his coordinate position and time. In general, the world function approach is useful in applications where 
high-accuracy measurements must be made over large distances. An application of recent interest is the design of 
space-based interferometers for precision sensing and surveillance purposes 15 lq| . In some of these designs, in order 
to achieve high-resolution imaging long base lines (hundreds of kilometers) must be used between Earth satellites and 
their separations may need to be accurate to within a micrometer or better. To unambiguously and accurately define 
such positions, a curved space-time approach should be used that takes into account the warping of the geometry of 
space-time due to gravitational effects. 

In section II, I write the generalization of Eq. (|l|) in terms of the world function and point out the limitations of 
navigation in curved space-time by electromagnetic beacons. In section III, I briefly describe the restricted problem of 
computing an observer's coordinate time if his or her spatial coordinates are known (a restricted type of navigation). 
In sections IV and V, I give the detailed equations applicable to navigation in the vicinity of the Earth in Schwarzschild 
coordinates and in coordinates that rotate with the Earth. 



II. NAVIGATION EQUATIONS 

The world function was initially introduced into tensor calculus by Ruse fl7|Jl8}] , Synge lisfl , Yano and Muto [po[ , 
and Schouten EM. It was further developed and extensively used by Synge in applications to problems dealing with 



measurement theory in general relativity |14|. In general, the world function has received little attention, so I give the 
following definition. Consider two points, Pi and P 2 , in a general space-time, connected by a unique geodesic path T 
given by a; 1 (it), where u\ < u < u 2 . A geodesic is defined by a class of special parameters u' that are related to one 
another by linear transformations u' = au + 6, where a and b are constants. Here, it is a particular parameter from 
the class of special parameters that define the geodesic T, and x l (u) satisfy the geodesic equations 

^ d ^ =Q (2) 



du 2 ^ k du du 

The world function between Pi and P2 is defined as the integral along T 

. 1 , . f" 2 dx l dxi , 

n(P 1 ,P 2 ) = -(u 2 -u 1 ) J 9iJ — —du (3) 

The value of the world function has a geometric meaning: it is one-half the square of the space-time distance between 
points Pi and P 2 . Its value depends only on the eight coordinates of the points Pi and P 2 . The value of the world 
function in Eq. (|^) is independent of the particular special parameter u in the sense that under a transformation from 
one special parameter u to another, u' , given by u = au' + 6, with x l {u) = x l (u(u')), the world function definition in 
Eq. (|J) has the same form (with u replaced by u'). 

The world function is a two-point invariant in the sense that it is invariant under independent transformation of 
coordinates at Pi and at P 2 . Consequently, the world function characterizes the space-time. For a given space-time, 
the world function between points Pi and P 2 has the same value independent of the metric-induced coordinates. A 
simple example of the world function is for Minkowski space-time, which is given by 

il(xi,4) = ^ri ij Ax i Ax j (4) 

where rjij is the Minkowski metric with only non-zero diagonal components (— 1,+1,+1,+1), and Ax 1 = (x l 2 — x\), 
2 = 0,1, 2, 3, where x\ and x\ are the coordinates of points Pi and P 2l respectively. Up to a sign, the world function 
gives one-half the square of the geometric measure in space-time. Calculations of the world function for specific 
space-times can be found in Refs. p4] , p2| -p4[ and application to Fermi coordinates in Synge (l4|] and Gambi et al. p5[ . 
The generalization to a curved space-time of the navigation Eq. (p]) is given by 



2 



n(xi,xi)=0, 5 = 1,2,3,4 



(5) 



where a;* = (i s , x s ) are the coordinates of the emission events, x l a = (t Q , x G ) are the observer coordinates and the world 
function is defined by Eq. (^|). Within the geometric optics approximation, Eq. (|5|) forms a natural generalization of 
the navigation Eq. (Q). In addition to Eq. (|5|), the appropriate causality conditions t s > t a or t a > t s , for s =1,. . .4, 
must be added. The set of relations in Eq. (^|) are manifestly covariant and invariant due to the transformation 
properties of the world function under independent space-time coordinate transformations at point Pi and at Pi. 

From the definition of the world function, the intrinsic limitations of navigation in a curved space-time are evident: 
the world function $l(P\,Pi) must be a single- valued function of P\ and Pi. In general, if two or more geodesies 
connect the points P\ and P2, then fi(Pi, Pi) will not be single- valued and the set of equations in Eq. may have 
multiple solutions or no solutions. Such conjugate points P\ and Pi are known to occur in applications to planetary 
orbits and in optics [Q. However, when the points P\ and P2 are close together in space and in time and the 
curvature of space-time is small, we expect the world function to be single valued and the solution of Eq. (||) to be 
unique. Therefore, navigation in curved space-time is limited by the possibility of determining a set of four unique 
null geodesies connecting four emission events to one reception event. In the case of strong gravitational fields as may 
exist in the vicinity of a black hole, or when the (satellite) radio beacons are at large distances from the observer in a 
space-time of small curvature, navigation by radio beacons may not be possible in principle. In such cases, one may 
have to supplement radio navigation by inertial techniques; see, for example, the discussion by Sedov f26|j . 



III. COORDINATE TIME AT A KNOWN SPATIAL POSITION 



Consider the restricted problem of an observer that knows his spatial position and wants to obtain his coordinate 
time 27 . A null geodesic connects the emission and reception events, so the value of the world function is zero, 

ft(t s ,x s ,i ,x ) = (6) 

where the satellite emission event coordinates (t s , x s ) and the observer spatial coordinates x D are known. Equation (|^) 
must be supplemented by the causality condition t Q > t s . The observer obtains his coordinate time by solving Eq. 
@fort . 

As a simple example of the application of Equation (^), consider an observer at a known spatial location who wants 
to compute his coordinate time in flat space-time in a rotating system of coordinates, y l , by receiving signals from a 
satellite at P s — (t s ,y s ). I take the transformation from Minkowski coordinates x l to rotating coordinates y % to be 
given by 



y Q = x° 



i fi\ 1 . n. o 

y = cos(— x ) x — sm(— x ) x 
c c 

y 2 = sin(— x°) x 1 + cos(—x°) x 2 

c c 

y 3 = x 3 (7) 

The world function is a two-point invariant that characterizes the space-time, so, in the rotating coordinates its value 
does not change. Using the world function for Minkowski space-time in Eq. (^) and the transformation to rotating 
coordinates y % in Eq. (pn, the world function is given by 

fi(P s ,P Q ) = f2(t s ,x s ,i ,x ) = 0(i s ,y s ,t o ,y o ) 
= \[{y -y s ?-c 2 {t -t s ) 2 ] 

+ {y s x y D ) • n sin(w(t - t s )) + 2 [y s ■ y D - (y s ■ n)(y Q • n)] sin 2 (-uj(t - t s )) (8) 

where n is a unit vector in the direction of the angular velocity vector to, which I take to be along the z-axis. I 
assume that the angular velocity of rotation is small, so that the time for light to travel from a satellite at y s to 
an observer at y G is small compared to the period of rotation 2-k/uo. I define the small dimensionless parameter 
5 = uj\y — y s |/c << 1. Equation (||) can then be solved for At — t Q — t s by iteration, leading to 

cAt = \y a - y s \ + -(y s x y ) • uj 

c 

1 " 1 ' [(ys X yo) ' + W 2 y s ■ y a - (y s ■ «)(y • w)] \ + 0(6*) (9) 



2c?"" " 3| y -y ( 
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The first term on the right side of Eq. (|9|) divided by c is the time for light to travel from the emission event at the 
satellite, (t s , Ys), to the observer at event, (t ,y ), in the absence of rotation. The second term is the celebrated 
Sagnac effect [£8-31 , which depends on the sense of rotation of the coordinates (sign of uS). The third term is a 
higher order correction that is independent of the sense of rotation; i.e., it is the same when w — > — u>. This term is 
on the order of 5xl0~ 14 s for satellite and Earth angular velocity of rotation parameters appropriate to the GPS. 
Equation leads to the standard expression for the Sagnac effect when we take the difference of propagation times 
for clockwise (AT-) and counterclockwise (AT+) propagation of light along the limit of a sequence of tangents on 
the perimeter of a circle, AT + — AT- = 4A ■ ui/c 2 , where A is the included area Note that the third term 

in Eq. ^ does not contribute to the difference of round-trip times, AT + — AT-, so it is not measurable in a Sagnac 
experiment. However, this term does contribute to a determination of coordinate time. 



IV. NAVIGATION IN THE VICINITY OF THE EARTH 



In the vicinity of the Earth, the gravitational potential can be approximated by |32j 



GM 



I ( ^ ) />2«-<.s(0)l 



(10) 



where p2(%) = (3a; 2 — l)/2 is the second Legendre polynomial and J2 is the Earth's quadrupole moment, whose value 
is approximately J2 = 1.0 x 10~ 3 . However, for navigation in the vicinity of the Earth, we can neglect J2 since it is 
three orders of magnitude smaller than the dimensionless coefficient of the monopole potential GM/r, which already 
contributes small corrections to propagation of electromagnetic radiation. I also neglect the effects of the rotation of 
the Earth, which give rise to small terms in the metric of space-time go a , since these effects are completely negligible 
at the present time |31|. Therefore, the Earth's gravitational field can be sufficiently accurately described using the 
Schwarzschild metric H33H 



ds 2 



1 - 



2GM \ 2 , 2 

c 2 dt 2 + 

c z r I 



dr 2 



1 - 



2GM 



l (de 2 



sin 



(11) 



In Eq. (y_l|) , I neglect the gravitational field of the sun and other planets, since the Earth is in free fall and these fields 
are essentially (up to tidal terms) cancelled as a result of the equivalence principle. 
Using the transformation to rectangular coordinates 



x u = ct 



r sin o cos < 
r sin 6 sin q 
r cos 9 



(12) 



and expanding in the small parameter GM/c 2 r, the metric for the Schwarzschild space-time can be written to first 
order as a sum of the Minkowski metric, r]ij, and the deviation from flatness tensor hij as 



where hij is given by 



ds 2 = gij dx l dx 3 = (rjij + hij) dx l dx 3 
(cdt) 2 (x a dx a f 



2GM 



(13) 



(14) 



The assumption that GM/c 2 r <C 1 is a restriction on the region of validity of Eq. ( |l3| ) to large r compared with the 
gravitational radius of the Earth, which is 2GM/c 2 ks 0.88cm. 

Following Synge, I approximate the world function for the metric in Eq. ( [T3| ) by replacing the integrals over the 
geodesic T by integrals along a straight line, and taking the special parameter u to vary in the range < u < 1, which 
leads to [Ol 



1 1 ■ f l 
Q(x\, x° 2 ) = - rjij Ax 1 Ax 3 + - Ax* Ax 3 / hijdu 

2 2 J 



(15) 
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I find that explicit evaluation of these integrals (see Appendix B) leads to p4 



x 2 - Xi 



GM . 

— 5- X2 - Xi (cos fell - COS 62 ) 



x 2 - xi 



log 



tan(^) 



(16) 



where cAt = — x{, and 6\ and 62 are defined by 



cosfl a = ^',? X2 Xl , } , a = 1,2 (17) 
|x a ||x 2 - Xl 

See Appendix C for an estimate of the error in Eq. (|ll|). The first term on the right side of Eq. ( |l6|) is the world 
function for Minkowski space-time, given in Eq. (^). The second and third terms in Eq. ( jlq ) give the corrections to 
the world function of Minkowski space-time due to the gravitational effects of mass M. The expression in Eq. Jl6| ) 
can be used in Eq. (||) as a basis for navigation, or in Eq. (|^) for computing coordinate time in the vicinity of the 
Earth. 

As an example of using the world function in Eq. (|l^), I consider determining the coordinate time of an observer at 
a known spatial position x G in the vicinity of Earth, using standard Schwarzschild coordinates. Taking the satellite 
position as x s = Xi, the observer position x = x 2 , and making use of the small parameter GM/(c 2 \x — x s \) <C 1, I 
solve Eq. (0) by iteration, leading to 



1 At = 



GM 



2 log- 




cos 0, 



cos ( 



(18) 



The first term in Eq. ( fL8| ) divided by c is the time for light to propagate from x s to x Q . The second term is the 
small correction due to the presence of the Earth's mass M distorting the space-time in its vicinity. This expression 
takes into account the delay of the electromagnetic signal in a gravitational field (see, for example, Ref. |35||Ip6|] and 
references cited therein). For a satellite directly overhead, where x D and x s are co-linear with the origin of coordinates, 
Eq. ( |l8| ) leads to the result 



cAt 



2GM , |x s 

— 2— 



(19) 



Equation Eq. ( |l9| ) is obtained by considering the spatial geometry shown in Figure ||. Define a triangle by the three 
points: origin at O, satellite at P s with coordinates x s , and observer at P with coordinates x Q . From point P a , draw 
a line P Q M perpendicular to OP s and define its length to be h. Equation ( JTg| ) is obtained from Eq. fllg| ) by setting 



6 S = 61, O = 02, x s = Xi, and x D = x 2 , using the definitions in Eq. ( |17| ) and the relations 

2 



cos 6 S 



cosc/„ = 



where 



-1 
-1 



h, 
1 



2/? 



h = (x s - X Q ) • — ^j- 



211 



(20) 
(21) 

(22) 

(23) 
(24) 



and taking the limit as h — > 0. Note that Eq. (18) and (|19|) are expressed in terms of Schwarzschild spatial coordinates 
(and do not contain the temporal coordinatesjof satellite x s and observer x , because the Schwarzschild space-time 
is static; i.e., the space-time admits a hypersurface orthogonal time- like Killing vector field |37t]. 
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V. NAVIGATION IN THE VICINITY OF THE EARTH IN ROTATING COORDINATES 



In practical navigation problems, an observer or user of a satellite navigation system is often interested in their 
space-time position with respect to the Earth-which defines a rotating coordinate system. To compute an observer 
space-time position in a rotating system of coordinates, I apply the navigation Eq. (||) in the rotating system. Having 
computed the world function for Schwarzschild space-time in Eq. ([l6]) in a 'nonrotating' system of coordinates x l , the 
invariant nature of the world function can be used to write the world function in a rotating system of coordinates, y l , 
using the transformation in Eq. (^): 



n{x\,x 3 2 ) = fi(*i,yi,* 2 ,y2) 

= \mj y l y J + af 

gm 1 
+- 



[(y 2 - yi ) 2 + 2AF] 

GM .. ., GM 

-— 2-(|yi| + |y2|) + — 2- 



— 2 [(y 2 - yi) 2 + (c At) 2 + 2 AF] log 



ft + A ) ^ ■ y* - AF i 

yi y2 ' 



where 



AF = (yi x y 2 ) • n sin(wAt) + 2 [yi • y 2 - (yi • n)(y 2 • n)] sin 2 

=yi-(y 2 -yi)-AF+| yi | [(y 2 - yi ) 2 + 2 AF] 1/2 

F 2 = y 2 • (y a - yi) + AF + |y 2 | [(y 2 - yi ) 2 + 2 AF] 1/2 



UjAt 



(25) 



(26) 

(27) 
(28) 



For navigation in the vicinity of the Earth in rotating coordinates, the observer must solve the four simultaneous 
Eqs. (o) using the world function in Eq. (|2^). In general, this must be done numerically. 

Consider the simpler problem of determining an observer's coordinate time at a known spatial location. Analytic 
results can be obtained in this case. I solve Eq. (||) for At = t Q — t s using the world function in Eq. (^) by defining 
three small dimensionless parameters, x = u>At, S = ui\y — Ys\/c, and a — GMuj/c 3 , and solving for i as a function 
of 6 by iteration. This leads to 



cAt = |y c - y s | + ~(y s x y ) • u> + -^\yo ~y s 



[(y« x y ) • ^] 2 
|y - y s | 2 



+ ^ 2 y s • y Q - (y« • w)(y D ■ w) 



GM 



2 , i (y» x y ) • u 
c |y -y s | 



log 



tan(^) 
tan(%). 



GM 



1 1 (ys x y p ) ■ uj 
c |y -y s | 



(cos 9 S — cos o ) 



(29) 



In Eq. (p9|), I have dropped small terms of 0(aS), 0(a 2 ), and 0(5 A ). Equation ( p9| ) gives the coordinate time for the 
signal to travel from the source to the observer, At = t a — t s . Since the propagation time At is given in a rotating 
system of coordinates, Eq. ( p9| ) contains the Sagnac effect, modified by the presence of mass M (the Earth). The first 
term corresponds to the propagation of the electromagnetic signal from the satellite at y s to observer at y Q in (flat, 
nonrotating) Minkowski space-time coordinates. The second term is the standard Sagnac correction term (expressed 
in terms of coordinate time) that depends on the sense of rotation uj, which appeared in rotating coordinates in vacuum 
(see Eq. ([)])). The third term on the right side of Eq. (^9|) also appears in flat space-time and, as previously mentioned, 
it does not depend on the sense of the rotation; it is the same when ui —> —uj and, hence, does not contribute to 
the Sagnac effect. The last two terms on the right-hand side of Eq. (|2^) are corrections to the coordinate time of 
propagation due to the Earth's mass M. Note that the coefficients of each of the last two terms also depend on uj. 
This represents a modification of the Sagnac effect due to the presence of mass M. 

The rotation of the coordinate system leads to a break in spherical three-dimensional symmetry. Note that terms 
two and three on the right side of Eq. (|29|) , which are due to the rotation, have a cylindrical symmetry determined by 
the direction of uj, as expected. On the other hand, the last two terms that depend on the mass M have coefficients 
that have a constant (spherically symmetric) term plus a term that depends on the sense of rotation (linear in uj). 
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VI. CONCLUSION 



If an observer simultaneously receives electromagnetic signals from four electromagnetic beacons in flat space-time, 
he or she can compute their position in space-time by solving the four light cone Eqs. (|l|). This procedure is routinely 
carried out everyday by users of satellite navigation systems such as GPS and GLONASS. Equation (jl]) neglects 
small effects of gravitational fields on electromagnetic signal propagation. In this work, I have included these effects 
in a natural way using the two-point invariant world function developed by J. L. Synge. I have given a simple 
covariant and invariant formulation of the n avig ation problem in Eq. (^|). An approximation to the world function 
in Schwarzschild coordinates is given in Eq. (16), and in coordinates that rotate with the Earth in Eq. (25). In the 
future, approximations to the world function may be obtained for the case of the Eddington and Clark metric |3S| ] 
and Eqs. (jlj) may serve as the basis for high-accuracy navigation throughout our solar system. 

I have implicitly used the geometric optics approximation by assuming that electromagnetic radiation propagates 
along null geodesies. Furthermore, when navigating inside the Earth's atmosphere, or when radiation traverses the 
atmosphere, there is an additional signal delay that is bigger than the effects discussed here and must be taken 
into account in detail. I have neglected these effects. Consequently, the results here are valid for observers in orbit 
around the Earth outside the Earth's atmosphere. Work is in progress to extend the world function approach to 
simultaneously include gravitational and index-of-refraction (atmospheric) effects. 

The world function approach used here can be applied to areas that require precise definitions of distances over 
satellite-to-ground length scales. A particular application area of recent interest, where high-accuracy position is 
required over hundreds of kilometers, is flying satellites in formation, for precise sensing and surveillance applications. 
The high-accuracies needed in definition of satellite positions may require a curved space-time theoretical framework 
that includes gravitational corrections to flat Minkowski space-time. 



APPENDIX A: CONVENTIONS AND NOTATION 



Where not explicitly stated otherwise, I use the convention that Roman indices, such as on space-time coordinates 
x 1 take the values i = 0, 1, 2, 3 and Greek indices take values a = 1, 2, 3. Summation is implied over the range of the 
index when the same index appears in a lower and upper position. In some cases, such as Eq. (|l4|), Greek indices 
summation is implied when both indices are in the upper position, such as in x a dx a . 

If x l and x % + dx l are two events along the world line of an ideal clock, then the square of the proper time 
interval between these events is dr — ds/c, where the measure ds is given in terms of the space-time metric as 
ds = —gijdx l dx :l . I choose gij to have the signature +2. When gij is diagonalized at any given space-time point, 
the elements can take the form of the Minkowski metric given by 7700 = — 1 ; Vafs = S a p, and r/oa = 0. 



APPENDIX B: INTEGRALS 



Two integrals are needed to explicitly evaluate the world function in Eq. ( |l5| ) (see Ref. [p4| ): 



■ du 



Xi - x 2 



[ tan(f ) 
log 1 



tan(* 



(Bl) 



Ax a A 



1 x a x fi 



du = Ixi - x 2 | 



log 



'tan(^)' 

.Ml). 



cos(6>i) — cos(6* 2 ) 



(B2) 



where 9\ and 62 are defined in Eq. @. 



APPENDIX C: ERROR IN APPROXIMATION OF THE WORLD FUNCTION 

In Eq. (|l6|), I use Synge's approximation to the world function for the metric in Eq. (|l3|), which entails replacing 
the integrals over the geodesic T by integrals along a straight line L 

1 f dx l dxi 1 f dx i dxi 

n(p u p 2 ) = -(u 2 - Ul ) gij ——du^-( U2 -u 1 ) g^^du (ci) 
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I take the straight line L to be given by 

x l L (u) — k{u2 — u)x\ + klu — Ui)x\ 



(C2) 



where u\ < u < u-i and k — (u,2 — wi) -1 . The error made in this approximation can be computed as follows. Consider 
two points, P\ and P%, connected by a family of curves x l (u, v), where u and v are independent parameters. I assume 
that Pi = x*(iti, v) and P2 — x l {u2, v), for all v. It is convenient to define the two vector fields 



IP 



dx l 
du ' 



V 



dx l 
dv 



(C3) 



Note that by construction V l (ui,v) — V l (u2,v) = 0. Furthermore, assume that a unique geodesic T connects the 
points Pi and P2, and that this geodesic T is given by the curve with a particular value of v = v , namely x l {u, v Q ) 
for ui < u < U2- Define the integral 



I(V) = ~(>2 -Ui) 



9ij 



U l U J du 



Now, assume that we have a space-time of small curvature and expand the integral I(v) about the geodesic T 



i(v) = n(Pi,p 2 ) 



dl(v ) 
dv 



(v -v ) + - 



1 d 2 I(v ) 



2 dv 2 



(V - V f + 



(C4) 



(C5) 



The second term, dl(v a ) / dv, vanishes since this is the definition of a geodesic curve. The error in replacing the integral 
over a geodesic by an integral over a nearby curve can then be estimated by the third term on the right side of Eq. 
(|C5|). Using the relations 



and 



I find the approximate error is given by 



Su Sv 



SU l (u,v ) 
Su 



(C6) 



(C7) 



i(v)-n(Pi,p 2 ) 



(V - V f 



1 d 2 I{v ) 

2 dv 2 

-(u 2 - ui)(v - v ) 2 



9ij 



SU' 1 5U J 
Sv Sv 



+ R abc dV a U b U c V c 



du 



where the Riemann tensor is given by 

R 

the afhne connection is 



jkm 



jm,k 



jk,m 



jm a k jk arn 



(C8) 



(C9) 



r )k = -jf (9jl,k + 9u,j - 9jk,i) 



(CIO) 



and ordinary partial derivatives with respect to the coordinates are indicated by commas. In Eq. (CS), I have dropped 
terms of order 0(v — v a ) 3 . 

To estimate the error incurred in Eq. by approximating the world function integral over a geodesic by an 
integral over the straight line in Eq. (C2), I construct an explicit parametrization of curves connecting Pi and P2: 



x 2 (u,v) = — — -x\(u) + —Xp(u) 
v„ v„ 



(Cll) 



where x l L (u) is given by Eq. ( |C2| ) and Xp(u) is a geodesic connecting the points Pi and P2. For simplicity, in Eq. 
( pll ) I am assuming that the parameter u = at Pi and u = U2 at P2- For the geodesic connecting the points Pi 
and P 2 , I use a series solution of Eq. (0) 
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(C12) 



where = x l 2 — x\. Using the curve parametrization in Eq. ( Cll ), and carrying out the required calc ulat ions, I find 
an estimate of the error in approximating the world function integral by a straight line, given by Eq. jC£p, to be 



%P a6 (p 1 )r^(p 1 )r^ d + o(3) 



(C13) 



where third order terms have been dropped, and T l ab (Pi) is the connection evaluated at Pi. I have taken ^) a b,c and 
lab,cd to be of O(l). 



APPENDIX D: NOTE ON ITERATIVE SOLUTION OF NAVIGATION EQUATIONS 



Equation (5|) is a set of four nonlinear algebraic equations for the four coordinates (t ,x ). A simple method of 
solution can be applied by linearizing and solving the system by iteration. Setting x 3 = x 3 (n) + 8x 3 {n) in Eq. (^) and 
expansion to first order in Sx 3 (n) gives a linear set of equations for 8x k (n) 



M^6x k (n) = -Q(xi,xi(n))=Q, 



a = 1,2,3,4 



(Dl) 



where 5x (n) is the correction to the n trial value x 3 (n) and 



ak 



an(xi,xi) 



dx 3 



-x J (n) 



(D2) 



I start the iteration by making an ansatz x 3 (l) for user coordinates x 3 , and solve Eqs. (Dl) for the first correction 
5x k (l). The improved solution is then taken to be x J Q (2) = x 3 (l) + Sx k (l) and substituted back into Eq. (Dl). 
Iteration is continued until sufficiently small corrections Sx k (n) are computed. 
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Figure 1 

FIG. 1. Space-time diagram showing two (of the four) satellites with world lines Si and S2. Electromagnetic signals are 
emitted at Pi and P% % and reach the observer at O. 
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Figure 2 

FIG. 2. The spatial geometry of satellite at Ps and receiver at Po, with coordinates x s and x , respectively, showing the 
angles 9 a and O . 
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